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Abstract 

Using a plane wave pseudopotential approach to density functional theory we investigate the 
electron localization length in various oxides. For this purpose, we first set up a theory of the band- 
by-band decomposition of this quantity, more complex than the decomposition of the spontaneous 
polarization (a related concept), because of the interband coupling. We show its interpretation 
in terms of Wannier functions and clarify the effect of the pseudopotential approximation. We 
treat the case of different oxides: BaO, a-PbO, BaTiC"3 and PbTi03. We also investigate the 
variation of the localization tensor during the ferroelectric phase transitions of BaTiC"3 as well as 
its relationship with the Born effective charges. 

PACS numbers: 
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I. INTRODUCTION 



In the study of periodic crystalline solids, the electronic ground-state wave function is 
usually described in terms of Bloch functions, delocalized on the whole system. As a conse- 
quence, for a long time, the understanding of electron localization in crystalline solids was 
mainly based on approximate pictures. Nevertheless, the basics of a quantitative character- 
ization of e.ectron locaWioa had a lrea d y oeen footed by W. KolmH in 1964. Receatl, 
this problematic was renewed, thanks to the development of the theory of polarization based 
on Berry phases!!®, and the rigorous definition of the position operator in periodic sys- 
tems!®!. These ideas have been further developped using a cumulant generating function 



approacnli 

Following these advances, Sgiarovello and co-workers0, have computed the localization 
lengths for different cubic semiconductors, in the framework of the Kohn-Sham density 
functional theory0'El (DFT). They showed that the degree of electron localization is quite 
different for the various investigated materials. These results encourage to pursue the study 
to other insulating crystals. 

The localization tensor (r a rp) c , can be computed from the periodic part of the Bloch 
functions u n k(r) and their first derivatives with respect to their wavevector : 

v. r , N 




(1) 
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where V c is the volume of the primitive unit cell, N the number of doubly occupied bands 
and a, (3 are two cartesian directions. The localization tensor is a global quantity that 
characterizes the occupied Kohn-Sham manifold as a whole (all k-points and all bands). 
This statement calls for two comments. First, applications of DFT to solids often make use 
of the frozen-core and pseudopotential approximations, while Eq. ([I]) requires an all-electron 
calculation. Second, the behavior of core and valence electrons is treated globally while both 
kinds of electrons are expected to exhibit strongly different localization properties interesting 
to identify independently. 

The localization tensor has been shown0 to give a lower bound for the spread of maxi- 
mally localized Wannier functions (WF) as defined by Marzari and VanderbiltBE! (hereafter 
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cited as MV) . In order to get some insight into the physics of the chemical bonds in molecules 
and solids, such WF are usually constructed considering only a restricted number of elec- 
tronic bands close to the Fermi level. The spread of the resulting WF is strongly dependent 
of the electronic states included in the minimization process. In this context, it seems inter- 
esting to try to identify the intrinsic localization of the electrons in a specific set of bands 
and to understand how this quantity is affected when including other bands. This would 
allow to solve the problem associated to the use of pseudopotentials and to characterize 
separately the behavior of core and valence electrons. 

This paper is organized as follows. In Sec. |l| we propose a decomposition of the lo- 
calization tensor into contributions originating from isolated sets of bands composing the 
energy spectrum of a solid. Using a simple model, we then illustrate the role of the cova- 
lent interactions on the different terms of the decomposition. We also make a connection 
between the localization tensor and the Born effective charges and we discuss the relation 



between pseudopotential and all-electron calculations. In Sec. I XT], we give the technical 
details underlying our first-principles calculations and we point out the differences between 
our method and that applied in Ref. |TT[ In Sec. ^ and M, we present the results obtained 
on two ferroelectric perovskites (BaTiOa and PbTiOs) as well as on two binary oxides (BaO 
and a-PbO). We investigate the variation of electron localization during the phase transi- 
tions of BaTiC>3 and show that the evolution is compatible with the electronic structure of 
this compound. 



II. BAND BY BAND DECOMPOSITION OF THE LOCALIZATION TENSOR 



A. Formalism 



Contrary to the polarization and the Born effective charges, for which band-by-band 
decompositions have been previously reported0'0S0, the localization tensor (Eq. (JTJ) ) 
involves scalar products between Bloch functions of different bands, making the identification 
of the contribution of isolated sets of bands less straightforward. In order to explain this 
fact, we have to remember that the localization tensor is related to the second moment of 
WF while the Born effective charges and the spontaneous polarization are linked to their 
first moment. From standard statistics, it is well known that these two quantities do not 
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add the same way : when considering two random variables X\ and x 2} the mean value of 
the sum x\ + X2 is simply the sum of the mean values while the variance of the sum is the 
sum of the variances plus an additional term, the covariance. 

These considerations can be transposed in the simple context of a confined model system 
made of two orthonormalized states ipi(x) and if^(x). The total many-body wavefunction 
ty(xi,x 2 ) is a Slater determinant constructed on the one-particle orbitals. The center of 
mass is given by the expectation value of the position operator X = J2i=i 2 ^ 

X=(*\X\*) = Y,y>i\&\ii>i) (2) 

i=l,2 

while the total spread (two times the localization tensor) is related to X 2 , 

a 2 = (^|X 2 |^) - (^|X|^) 2 

= ^[(V# 2 |^) - -2{^ 1 \x\^ 2 )^ 2 \x\^ 1 ). (3) 

i=l,2 

We see that the first moments of the one-particle orbitals add to form the total dipole of 
the many-body wavefunction. On the contrary, the total spread is not equal to the sum of 
the individual spreads of ip\ and ip 2 but involves also matrix elements of the one-particle 
position operator x between tpi and i\) 2 . The additional term would be absent if the many- 
body wavefunction was a simple product of the one-particle orbitals. It arises from the 
anti-symmetry requirement. In analogy with the language of statistics, we will name it the 
covariance. 

Based on these arguments, we can now define a band-by-band decomposition of Eq. (p. 
Suppose that the band structure is formed of N g groups labelled Gi, each of them composed 
of rij bands (i = 1, N g ). The variance of a particular group Gi is defined as 



(r a rp) c {Gi) = 




where the sums have to be taken over the bands of group Gi. The covariance of two groups 
Gi and Gj (i 7^ j) is given by the following relationship: 



(r a rp) e (Gi,Gj) = 
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Using these definitions, the total tensor, associated to the whole set of occupied bands, can 
be written as 

(r a rp) c = Jj^2 n i\ ( r a r fi)c(Gi) + ^ j n j {r ol r l3 ) c (Gi 1 G j ) \ . (6) 

i=l I j+i J 

The variance (r a r^) c (Gj) is intrinsic to an isolated set of bands. It is related00 to the 
quantity Qj introduced by MV through the relation 

3 

&i = rii 22(r a r a ) c (Gi). (7) 

a=l 

In a one-dimensional crystal, Qj is simply the lower bound of the total spread Q of the 
WF§, 

0= X> 2 >„-<r£], (8) 

that can be realized by choosing an adequate phase factor for the Bloch functions. In a 
three-dimensional crystal, it is no more possible to construct WF that are simultaneously 
maximally localized in all cartesian directions. It is only possible to minimize their spread 
in one given direction as realized for the so-called hermaphrodite orbitals introduced in 
Ref. [TT] : these particular functions are localized (Wannier-like) in a given direction a and 
delocalized (Bloch-like) in the two others. The variance of a particular group of bands 
(r a r a ) c (Gi) is the lower bound of the average spread ^- J2 n eGi l( r a)n ~ ( r a)«] where the sum 
is taken over all Wannier-like functions in the unit cell belonging to group Gj. This lower 
bound is reached for WF that are maximally localized in direction a. The variance therefore 
gives some insight on the localization of the electrons within a specific set of bands taken 
independently. This localization is affected by the hybridizations between atomic orbitals 
giving rise to the formation of the considered electronic bands within the solid so that the 
variance can act as a probe to characterize these hybridizations. 

The covariance is no more related to an isolated set of bands. It teaches us how the 
construction of WF including other bands can improve the localization. As discussed by 
MV, the definition of groups of bands in a solid is not unique and sometimes there is a 
doubt about which bands have to be considered together. If we consider two sets of bands 
Gi and Gj as one single group, its total variance is the sum of the individual variances and 
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covariances, that have to be rescaled by the number of bands in each group 

(r a rp) c = — - — {rii [(r a r p ) C (G + nj{r a rp) c (Gi, Gj)\ +rij [(r a rp) c {Gj) + ni(r a rp) c (G j ,G i ))} . 

(9) 

Until now, we considered separately the two cartesian directions a and (3. Stronger 
results can be obtained when diagonal elements of the localisation tensor are considered, or 
when this localication tensor is diagonalized, and the eigenvalues are considered. Different 
inequalities can be derived. In particular, from Eq. (|^), it appears that the covariances for 
a — (5 are always negative. This means that the diagonal elements of the full tensor are 
always smaller than those obtained by the sum of the diagonal variances. In other words, it 
is always possible to obtain more strongly localized orbitals by constructing WF considering 
more than one group of bands. As a consequence the covariance appears as a tool to identify 
which bands have to be considered together in the construction of WF in order to improve 
their localization. 

In appendix [A] we give an interpretation of the variance and covariance in terms of the 
optical conductivity. It illustrates from a different viewpoint the influence of the fermionic 
nature of the electrons on the localization tensor : the appearance of the covariance in Eq. 
(|) is a direct consequence of the Pauli principle. 

B. Simple model 

In this section we will investigate a one-dimensional model system. This will help us 
to understand the role of the covalent interactions on the electron localization length and 
related quantities like the Born effective charges. We will deal with a confined system for 
which the localization tensor can be computed from matrix elements of the position operator 
and its square as described in Ref. [□]. 

Let us consider a diatomic molecule XY. In order to describe the chemical bonds of this 
model system we adopt a tight-binding schemel! defined by the hopping integral t and the 
on-site terms A and —A. We will call a the interatomic distance and ipx, ipY the s-like 
atomic orbitals that are used as basis functions. The hamiltonian can be rescaled by A 
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(A=t/A) in order to become a one parameter hamiltonian defined by 

*=(-;;)■ ™ 

We further assume that ipx is centered at the origin, ipy in a and that these two functions 
do not overlap at any x 

tpx{x)Mx-a) = 0. (11) 
The eigenfunctions of the hamiltonian correspond to 

0i i2 (x) = u lt2 ipx(x) + v lt2 ipY(x - a) (12) 
where the coefficients Ui^ and t> lj2 can be expressed in terms of the bond polarity^! 



/ 1— a 



. l—a p / l+a p 



(13) 



In order to see the meaning of the different terms appearing in the band by band de- 
composition of the localization tensor and the Born effective charges let us first consider the 
molecular orbitals independently. 

The variance of state <f>% can be computed from the coefficients ui and v±. It writes 

<* 2 > c (i) = 4^ + 4^ + (i4) 

where a\ and a\ are the second central moments of ipx and ipy ■ The variance of (f>2 is given 
by a similar expression. This quantity is composed of three positive terms that summarize 
the mechanisms that are able to delocalize the electrons with respect to the atomic orbitals. 
On one hand, the electronic cloud on a particular atom is not a delta-Dirac function but 
presents a degree of derealization related to a\ and a Y (first and second term). When 
the state (f>i is made entirely of ipx, that is, when a p equals one, the localization length is 
correctly equal to a x (first term). Incorporating more ipy changes the localization length 
in proportion of a p (the balance between first and second terms). On the other hand, the 
electrons can occupy two sites X and Y that are separated by a distance a (third term). 
This term scales as a 2 . Even a small covalent interaction is thus able to induce an important 
derealization if it acts on a sufficiently large distance. 



The Born effective charge of atom X is defined as the derivative of the dipole moment 
p with respect to a. This dipole moment is the sum of the nuclear and static electronic 
charges multiplied by the interatomic distance. The contribution coming from the electrons 
occupying state 0i is equal to 

Pi = —2eu\a = — e(l + a p )a (15) 



X,l = -fo= ~< l + ^ + ea (l + ^)3/ 2 ^ (16) 



where e is the module of the electronic charge. The derivative of Eq. (15) with respect to a 
gives the contribution of these electrons to the total effective charge 

7 dPi „ A OA 

Zj 

The first term is the (static) effective atomic charge! of atom X while the second term 

represents an additional dynamical contribution due to a transfer on electrons between X 

and Y during a relative atomic displacement. The contribution of the electrons occupying 

state 02 is given by a similar expression 

* dp 2 A dA 

Z ^ = to = - 6(1 - Qp) ~ ea (l + A^/^ 

This simple model illustrates how both the variance of the localization tensor and the 
Born effective charges depend on the covalent interactions defined by the parameter A. 
The variance is a static quantity depending on the amplitude of the covalent interations 
only while the the Born effective charges are dynamical quantities that also depend on the 
variations of these interactions during a relative atomic displacement. 

If we now consider the states 0i and 02 as a single group we have to add their variances 
and covariances to get the whole localization tensor. The covariance reduces to 



By adding this covariance to the variance in Eq. fll4|) , we remove in some sense the de- 
localization induced by the covalent interactions. The total localization tensor becomes 
independent of the hopping A and the interatomic distance a. It reduces to the mean 
spread of the atomic orbitals ijjx and ipy '■ 

<* 2 > c = ^y^- (19) 

Eq. (|19D defines the mean spread of the WF constructed as linear combinations of 0i and 02 
that minimize the spread functional Q (see Eq. (H)). As shown by MV, these WF diagonalize 
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the position operator x projected on the subspace of occupied states. They are thus equal to 
the atomic orbitals since the hypothesis of zero overlap (Eq. (|TTD) implies (ipx\x\ipY) = 0. 

The total Born effective charge of atom X can be obtained by adding the nuclear charge 
Z* ore = 2e to the terms ( |TE| ) and flT7|). It is easy to check that for this model Z* x is equal to 
zero. This result can be interpreted in two ways. The point of view usually adopted is to 
say that the two molecular orbitals are of the opposite polarity so that the total dipole of 
the molecule vanishes. Based on the results of the preceeding paragraph, we can also affirm 
that each maximally localized WF is confined on a single atom so that no interatomic charge 
transfer can take place. 

This result suggests that the variance gives more informations about the localization of 
electrons of particular chemical bonds than the total localization tensor. It also illustrates 
the observation of Ghosez et that anomalous effective charges mainly come from 

hybridizations between occupied and unoccupied states. In fact, the different chemical 
bonds generate opposite effects so that a net charge transfer is possible only if some of them 
are unoccupied. 

In summary, we have illustrated the mechanisms that govern the variance of the local- 
ization tensor and the Born effective charges in the particular case of a one dimensional 
model system. The observations made in this section give us an intuitive understanding of 
how delocalized electrons can generate anomalous effective charges. Hybridizations between 
occupied states generate opposite effects that tend to cancel out when they are summed. 
Because of the simplicity of the above adopted picture, we have however to be careful when 
we apply this model to real materials. First, we considered only hybridizations between 
two types of atomic orbitals, while the chemical bonds in real systems generally result from 
more complicated interactions. In particular, we neglected on-site hybridizations that are 
also able to generate anomalous effective charges but that induce a stronger localization on 
the electronic cloud. Second, the hypothesis of zero overlap fllTD is not always fullfilled so 
that maximally localized WF constructed on the whole set of occupied states generally not 
reduce to the atomic orbitals. Nevertheless, this simple model will allow us to interpret 
some results in Sections [TV] and |V|. 
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C. Pseudopotentials 



As mentioned in the introduction, there is a fundamental problem in the computation of 
the total localization tensor when pseudopotentials are used. This is due to the fact that 
the localization tensor is related to the bands of the system as a whole : first, there is no 
cancellation between the core electrons and the nuclear charge, as it is the case in the com- 
putation of the total polarization; second, the localization tensor is a kind of mean over all 
bands, that combines strongly localized (core) states, and weakly localized (valence) states. 
This is clearly seen in Eq. ([!]), where the number of bands explicitely appears both as the 
denominator of the prefactor and in the two summations. The band-by-band decomposition 
allows us to overcome this problem partly, by focusing only on the variances of isolated 
groups of bands. Thanks to Eq. @ it is also possible to get some insight into the physics of 
the all-electron localization tensor when pseudopotentials are used. In this section, we focus 
on the diagonal elements of the electron localization tensor a = (3 (of course, any direction 
can be chosen as a). 

In an all-electron calculation, let us consider separately two sets of bands: core bands 
(labelled as 'co'), and valence bands (labelled as 'va'). The total localization tensor can be 
obtained from the localization tensors of each group of bands, combined with the covariance 
between the two groups of bands: 

(r a r a ) e = {n co (r a r a ) c (co) + n va (r a r a ) c (va) + 2n co n va (r a r a ) c (co,va)} . (20) 

Both variances (r a r a ) c (co) and (r a r a ) c (va) are positive quantities. The covariance times 
the product of the number of bands n co n va (r a r a ) c (co, va) , a negative quantity, must always 
be smaller in magnitude than each of the related variances multiplied by the corresponding 
number of bands. This translates to bounds on the diagonal elements of the total localization 
tensor : 

\n va {r a r a ) c (va) - n co (r a r a ) c (co)\ . n va (r a r a ) c (va) + n co (r a r a ) c (co) 

— \ r a r a)c _ • l^-U 

Tb co "T Tl va T^co < ""ua 

In the frozen-core approximation, (r a r a ) c (co) can be obtained from separate all-electron 
calculations for each atom of the system. The localization tensor of the valence bands is 
(likely) computed accurately in the pseudopotential approximation : the spread of the Wan- 
nier functions should be quite similar if estimated from all-electron valence wavefunctions 
or from pseudo- wavefunctions. 

10 



Thus, a bound on the diagonal elements of the localization tensor can be computed from the 
atomic wavefunctions of the core electrons and the pseudo-valence wavef unctions. In order 
to compute the covariance more accurately it is necessary to reconstruct the all-electron 
wavefunctions. This could be done following the ideas exposed in Ref. 133. 



III. METHOD AND IMPLEMENTATION 



In the remaining part of this paper, we apply the previous formalism to various ox- 
ides. The electronic wavefunctions are obtained within and the local density ap- 
proximation (LDA) thanks to the ABINIT0 package. At variance with a previous work on 
semiconductorsS, the first derivatives of the wave functions with respect to their wavevec- 
tor are not computed from finite differences but from a linear-response approach@ within 
the parallel-transport gauge. The wave functions are further transformed to the diagonal 
Kaug Both ground-state and first-order wavefunctions are expanded in plane waves up 
to a kinetic-energy cutoff of 45 Hartree. Integrations over the BZ are replaced by sums 
over a 8 x 8 x 8 mesh of special k-pointsi3. With these parameters, the convergence of the 
localization tensor for the investigated compounds is better than 10~ 3 Bohr 2 . The ionic- 
core electron potentials of the Ba, Pb, Ti and O atoms are replaced by ab initio, separable, 
extended norm-conserving pseudopotentials, as proposed by M. Teterll. Ba 5s, 5p and 6s 
electrons, Pb 6s, 5d and 6p electrons, Ti 3s, 3p and 3d electrons, O 2s and 2p electrons are 
considered as valence states. Beside calculating the localization tensor on bulk-materials, 
we also computed it on the isolated atomic systems Ba 2+ , Pb 2+ and O by placing each atom 
at the origin of a periodic supercell of 20 Bohrs. 

As shown by Sgiarovello et a/.0, the localization tensor and thus the variances and co- 
variances, are real. Moreover, they are obviously symmetric in a and (3. Consequently there 
exists a set of cartesian axes where they are diagonal and their eigenvalues are also real 
numbers. In the discussion of our results we will always work in this particular frame so 
that we do not need to consider the off-diagonal elements of the localization tensor. 
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IV. RESULTS 



A. Structural and electronic properties 

We will consider the two binary oxides BaO and a-PbO as well as the ferroelectric per- 
ovskites BaTiOs and PbTiOs. BaO has a rocksalt structure while the tetragonal a phase 
of lead oxide is formed of parallel layers of Pb and O atoms. BaTiOs and PbTi03 have 
a high-temperature cubic perovskite structure with five atoms per unit cell. As the tem- 
perature is lowered, the former compound undergoes a sequence of three ferroelectric phase 
transitions transforming to tetragonal, orthorhombic and rhombohedral structures while the 
latter compound only undergoes one single transition from the cubic to the tetragonal phase. 
We will consider explicitely the cubic, tetragonal and rhombohedral phases of BaTiOs as 
well as the cubic phase of PbTi0 3 . 

The electronic structures of these compounds have been previously studiedBSiE and 
are illustrated in Fig. p]. They are formed of well separated groups of bands. Each of them 
has a marked dominant orbital character and can be labeled by the name of the atomic orbital 
that mainly composes the energy state in the solid. The bands at the Fermi level are mainly 
composed of O 2p states that show significant interactions with other atomic orbitals like 
the well known O 2p-Ti 3d hybridization in BaTiOs an d PbTi0 3 . The bandstructures in the 
ferroelectric phases of BaTiOs are similar to that of the cubic phase. The phase transitions 
principally affect the bandgap and the spread of the O 2p bands while the positions of the 
deeper lying bands remain quite constant. The main difference in the electronic structures 
of BaO and BaTiOs on one hand and PbO and PbTiOs on the other hand comes from the 
presence or absence of Pb 6s electrons (that form the so called lone-pair in PbO). These 
electrons show a strong hybridization with the O 2p states. As a consequence the O 2p 
and Pb 6s bands are degenerate at the R point in PbTi0 3 and around the Z point in PbO. 
Consequently, we have to consider them as one single group of bands in the decomposition 
of the localization tensor. 

B. Localization tensor and Born effective charges 

As the total localization tensor is meaningless in pseudopotential calculations that do not 
include covariances with the core states, we focus on the variances of the different groups 
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FIG. 1: Band structures of BaO, cubic BaTi03, cubic PbTiOs and a-PbO. 

of bands. The values can be found in the Tables | and |I| where they are compared to the 
variances of the dominant atomic orbitals. We do not report any values associated to the 
deepest lying Ti 3s and Ti 3p bands although they have been included in our pseudopotential 
calculation. Their variances are in fact close to the atomic ones and they do not show any 
sizeable covariance with other bands in both BaTiC>3 and PbTiC>3. 

TABLE I: Variances (Bohr 2 ) of the Ba 5s, O 2s, Ba 5p and O 2p bands for the isolated atomic 
systems Ba 2+ and O, BaO and the cubic (C), tetragonal (T) and rhombohedral (R) phases of 
BaTi0 3 . 

System Str. Element Band 







Ba 5s 


O 2s 


Ba 5p 


O 2p 


Atom — 


(r 2 )c 


1.011 


0.929 


1.370 




BaO 


(r 2 )c 


1.065 


1.552 


2.023 


2.199 


BaTi0 3 C 


(r 2 )c 


1.091 


0.950 


2.189 


1.875 


T 


(rlh 


1.091 


0.945 


2.180 


1.852 




(rpc 


1.088 


0.965 


2.175 


1.842 


R 


(rl)c 


1.092 


0.963 


2.196 


1.862 




(rpc 


1.092 


0.949 


2.189 


1.804 



In the cubic crystals BaO, BaTiC>3 and PbTiC>3 as well as in the atomic systems, the 
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TABLE II: Variances (Bohr 2 ) of the 2s, Pb 5d and Pb 6s +0 2p bands in PbTi0 3 , a-PbO and 
for the isolated atomic systems Pb 2+ and O. 



System 


Element 


2s 


Band 
Pb 5d 


Pb 6s + O 2p 


Atom 


(r 2 )c 


0.929 


0.657 




PbTi0 3 


(r 2 )c 


1.874 


1.490 


1.749 


PbO 


(rl)c 


2.234 


1.142 


2.178 




(rpc 


1.724 


0.990 


1.968 



reported tensors are isotropic so that we only mention their principal values (r 2 ) c - This 
is no more true in the ferroelectric phases of BaTiC>3 where a weak anisotropy can be ob- 
served. The tensors have an uniaxial character as the corresponding dielectric ones: they 
are diagonal when expressed in the principal axes and the elements (r\) c and (r 2 ) c refer to 
cartesian directions perpendicular and parallel to the optical axis (that has the direction of 
the spontaneous polarization). A much stronger anisotropy is observed in a-PbO where the 
localization tensor has the same symmetry as in the ferroelectric phases of BaTiC>3. Due to 
its particular structure formed of atomic Pb-0 planes the electrons of each group of bands 
are more delocalized in a direction parallel ((rj_) c ) to the atomic planesS than perpendic- 
ular {(r 2 ) c ) to them. This observation agrees with our intuitive picture that the covalent 
interactions between atoms inside a layer are stronger than between atoms of different layers. 

Examining the variances of the different groups of bands we see that the Ba 5s electrons 
show a similar degree of localization both in BaO and BaTiC>3 also equivalent to that of the 
corresponding atomic orbital. On the contrary, the O 2s electrons behave differently in the 
materials under investigation: in BaTiC>3, their variance is close to the atomic one while they 
show a significant larger derealization in the three other compounds. It is in fact surprising 
to see the degree of derealization of the inner bands like the O 2s, Ba 5p or Pb 5d bands. 
In some cases like BaTiC>3, the electrons of these bands are even more strongly delocalized 
than those of the bands at the Fermi level. These results suggest that the corresponding 
atomic orbitals are chemically not inert but present non negligible covalent interactions. 
An interesting observation can be made for the O 2s and Pb 5d bands in PbTiC>3 and 
a-PbO. The derealization induced by the covalent interactions that generate these bands 
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tends to disappear when we consider them as one single group. In order to compute the 
variance of the whole O 2s and Pb 5d bands, we have to use Eq. @. As an example let us 
consider PbTiC>3. The different elements can be summarized in a matrix where the diagonal 
elements are the variances (Bohr 2 ) and the off-diagonal elements the covariances (Bohr 2 ) of 
the individual groups 

1.874 -0.240 
-0.240 1.490 

The total variance of the (O 2s + Pb 5d) group considered as a whole reduces to 0.734 Bohr 2 . 
For a-PbO, we obtain similar values of 0.732 Bohr 2 for (r\) c and 0.701 Bohr 2 for (r 2 ) c . These 
values can be compared to the mean spread of the atomic orbitals | (0.929 + 5 x 0.657) = 
0.702 Bohr 2 . 

The results presented above suggest that inner orbitals like O 2s, Ba 5p or Pb 5d are 
chemically not inert in the materials under investigation. This observation seems in con- 
tradiction with the conclusions drawn from partial density of states analysis^ that these 
states are rather inert. Nevertheless the inspection of the Born effective charges in BaO or 
BaTiO,0H confirms our observations that will now be illustrated for a-PbO and PbTiCh. 
This points out that the global shape of the bandstructure is less sensitive to the underly- 
ing covalent interactions than the variance of the localization tensor or the Born effective 
charges. 

In order to investigate the connection between the localization tensor and the Born effec- 



tive charges we report in Table |IT| the band by band decomposition of Z* Pb in PbTiC>3 and 
a-PbO. In the perovskite, this tensor is isotropic while in a-PbO it has the same symmetry 
as the localization tensor. The contribution of each group of bands has been separated into 
a reference nominal value and an anomalous chargefi For a-PbO, we observe the same 
anisotropy as for the localization tensor: the covalent interactions inside an atomic layer 
(Zp bJ _) generate larger anomalous contributions than the interactions involving atoms of 
different layers (Z Pb ^). By looking at the O 2s and Pb 5d bands we see that they gener- 
ate important anomalous charges that confirm our observations concerning the variances 
of these bands. Interestingly, in both materials these contributions cancel out when they 
are summed. We observe thus the same tendencies for the Born effective charges and the 
localization tensor: the effects induced by the covalent interactions between inner orbitals 
tend to disappear when the resulting bands are considered together. 
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TABLE III: Band by band decomposition of the Born effective charges (a. u. of charge) in 
PbTiOa and a-PbO. The contributions have been separated into a reference nominal value and an 
anomalous charge. 





PbTiOs 


a-PbO 




Band 




^pb± 


*Pb\\ 


Core 


14.00 


14.00 


14.00 


2s 


+ 3.47 


+ 1.89 


+ 0.26 


Pb 5d 


-10 - 3.36 


-10 - 1.80 


-10 - 0.40 


Pb 6s + 2p 


-2 + 1.78 


-2 + 1.06 


-2 + 0.48 


Tot. 


2 + 1.89 


2 + 1.15 


2 + 0.34 



V. DISCUSSIONS 



Based on the simple model exposed in Sec. [11 B| we can suggest the following mechanism 
to explain the results presented in the preceeding section. The atomic orbitals O 2s and 
Pb 5d (for which the hypothesis of zero overlap (|TT]) is reasonable) present weak covalent 
interactions that generate the corresponding energy bands in PbTiOs and cv-PbO. When 
we construct maximally localized WF for each individual group, the resulting orbitals are 
delocalized on Pb and O atoms so that during an atomic displacement an interatomic transfer 
of charges — generating anomalous Born effective charges — is possible. The fact that the 
variance of the global (O 2s + Pb 5d) group of bands is close to the mean spread of the 
atomic orbitals suggests that the maximally localized WF constructed on these bands are 
similar to the original atomic orbitals. In other words, they are confined on a single atom. 
This confinement also suppresses the interatomic charge transfer so that the anomalous 
charges disappear. We can make similar observations for the Ba 5p and O 2s bands in BaO 
and BaTiOs, although, in the latter compound, the cancellation in the Born effective charges 
and the variance is not as complete as in the three remaining ones. This suggests that in 
the lead oxides as well as in BaO, the inner bands Pb 5d and O 2s (resp. Ba 5p and O 2s) 
mainly result from hybridizations between two types of atomic orbitals. At the opposite, in 
BaTiOs the Ba 5p and O 2s bands are formed of more than two types of atomic orbitals. 

Looking now at the bands at the Fermi level, we see that their variance is significantly 
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larger in BaO and a-PbO than in the corresponding perovskites and that it remains nearly 
constant in the different phases of BaTiOs. This latter observation seems surprising for two 
reasons, (i) The LDA bandgap presents drastic changes when passing from the cubic (1.72 
eV) to the rhombohedral (2.29 eV) phase. This increase suggests a much stronger localization 
of the O 2p electrons in the ferroelectric phases, (ii) The giant Born effective charges observed 
in the paraelectric phase0'0 imply an important reorganization of the electronic cloud during 
an atomic displacement. It appears surprising that this reorganization has such small effects 
on the localization tensor. These small variations are not restricted to BaTiOs but similar 
observations have been made in other ferroelectric compounds like LiNbO30. 

Considering point (i), we note that the correlation between the bandgap and the localiza- 
tion tensor is not as tight as one might think. The variance of the O 2p bands for instance is 
significantly larger in BaO than in BaTi03 in spite of the fact that its LDA bandgap (1.69 
eV) is close to the gap in the cubic phase of BaTi0 3 . 

Considering point (ii), we note that it is possible to have an important reorganization 
of the electronic charge without affecting the localization tensor a lot. Following the ideas 
of the Harrison modelil, the giant effective charges in perovskite ferroelectrics result from 
dynamical orbital hybridizations changes generating interatomic transfers of charges. In Fig. 
(a) we have drawn shematically an O centered WF in the cubic phase of BaTi0 3 along 
a Ti - O chain. Due to the O 2p - Ti 3d hybridization, this WF has a finite probability 
on the neighbouring Tii and Ti2 atoms. According to the Harrison model, a fraction of 




dt 



FIG. 2: Oxygen centered WF in the cubic phase (solid line) of BaTiOs ( a ) an d its variation during 
the transition to the tetragonal phase (dashed line) (b). 
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electrons is transferred from Tii to Ti2 during a displacement dr of the O atom (Fig. 0(b)). 
Even if the quantity of charges involved in this process is small, the large scale on which 
this transfer takes place (of the order of the lattice parameter) implies a shift of the WF 
center larger than the underlying atomic displacement and explains the anomalous effective 
charges. During the transition from the cubic to the tetragonal phase, the central O atom 
is displaced by few percents of the lattice constant a (^ = 0.045) with respect to Tii and 
Ti 2 . The resulting shift of the WF center generates the spontaneous polarization in the 
ferroelectric phase. 

Based on this simple picture the origin of the small variations of the O 2p variance during 
the phase transitions becomes more obvious: When the electrons are transferred from Tix 
to Ti 2 their distance to the initial WF center remains unaffected and their distance to 
the displaced WF center slightly decreases due to its shift towards Ti 2 . Mathematically 
speaking, due to the fact that the variations do not depend on the direction of the atomic 
displacement, they are of the second order in — . 

In order to get a numerical estimate of the charges transferred during this process and 
its impact on the localization tensor we can consider a one dimensional model WF whose 
square is the sum of three delta-Dirac functions 



This model only takes into account the derealization of the electrons on different atoms 
(third term of Eq. (|i~4[) ) while it completely neglects the derealization of the electronic 
cloud on the individual atoms (first and second term of Eq. (0)). In this particular case 
we can identify the localization tensor to the second moment of the WF. This is no more 
true in a real, three dimensional crystal. In BaTiOs for instance, the O 2p group contains 9 
different WF per unit cell located on three different O atoms. These WF extend in different 
spatial directions so that their average spread in the x-direction is lower than the spread of 
one single WF as the one shown in Fig. [| 

In Eq. (|22|), Z' Q represents the probability of the electrons to be found on the O atom. It 
can be computed from the value of the O 2p variance in the paraelectric phase of BaTiC>3 and 
the lattice constant a using the relation J x 2 \W n (x)\ 2 dx = (r 2 ) Cj o2 P - This yields Z' = 1.73. 
This quantity allows an estimate of the static charge of the O atom in BaTiC>3 by substracting 
three times Z' Q from the charge due to the nucleus and the core electrons O Is and O 2s. 




(22) 
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This yields Z , st = 4 - 3 • 1.73 = -1.19 e. 

When the O atom is displaced, the shift of the WF center is directly related to the 
quantity of charges e transferred from Tix to Ti 2 . The value of e can be computed from 
the value of the effective charge generated by the O 2p electrons {Zq 2 = —9.31) in the 
cubic phase0 by taking into account that the anomalous charges are generated by three WF 
located on the same O atom0. To get the polarization due to one single WF, we have to 
divide this quantity by 3 since each of them brings a similar contribution to Zq 2 . In the 
tetragonal phase, the model WF writes 

\W n (x)\ 2 = \ | 2 ~^~^ (x + a) + Z' S(x - dr) + 2 ~ Z ° + £ 5{x - a)} . (23) 

By identifying twice its first moment to ZQ 2p dr/3 one gets e = 0.0614 at the transition from 
the cubic to the tetragonal state. It implies a decrease in the spread of the model WF of 
0.18 Bohr 2 . 

This variation is larger than the observed one (0.023 Bohr 2 ). Part of the discrepancy is 
probably due to the fact that we considered Z* 02v to be constant along the path of atomic 
displacement from the paraelectric to the ferroelectric phase. Using the value of Zq 2 in 
the tetragonal phase we obtain a value of 0.0467 for e while the variance decreases of 0.12 
Bohr 2 . Moreover, one has to bear in mind that the localization tensor in BaTiC>3 is an 
average value that has to be taken over 9 WF. Six of them are centered on O atoms that ly 
in a plane perpendicular to the direction of the spontaneous polarization. They are probably 
less affected by the phase transition. As a consequence, the variation of the WF located 
on the remaining O atom (the one represented on Fig. |2|) is expected to be larger than the 
variation of the localization tensor. 

In summary, even if there is no formal connection between the real WF in BaTiC>3 and Eq. 
(p2|), this simple model shows that small variations of the localization tensor are compatible 
with giant effective charges and their interpretation in terms of the Harrison model. As 
illustrated with the model WF, the transfer of charges along the Ti— O chains only implies 
a slight decrease in the spread of one single WF. This decrease is expected to be larger than 
the decrease in the variance because this latter quantity is an average value over 9 WF that 
are not modified to the same extent during the phase transition. 
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VI. CONCLUSIONS 



Using a plane wave-pseudopotential approach to DFT we computed the electron local- 
ization tensor for various oxides. Our study was based on the work on semiconductors 
performed by Sgiarovello and co-workers but used linear-response techniques to compute 
the first-order wavef unctions. 

In order to investigate the properties of electrons occupying individual groups of bands 
independently, we first set-up a band by band decomposition of the localization tensor. In 
analogy with the field of statistics we had to distinguish between variance and covariance in 
this decompositon. The significance of these new concepts was illustrated in terms of WF and 
explained on a simple model. The variance allows to get some insight into the hybridizations 
of atomic orbitals. The covariance can be useful to help constructing maximally localized 
WF: It identifies the bands that have to be considered together in order to improve their 
localization. We also made a connection between the localization tensor and the Born 
effective charges and we discussed the difference between all-electron and pseudopotential 
calculations. 

We applied these techniques to binary oxides (BaO and a-PbO) as well as perovskite 
ferroelectrics (BaTi0 3 and PbTi0 3 ). By considering first the electrons of the inner bands 
we showed that some of them present a strong derealization with respect to the situation 
in an isolated atom. This observation suggests that the underlying atomic orbitals are 
chemically not inert but present non negligible covalent interactions. This fact had been 
confirmed from an inspection of the Born effective charges. 

Finally, the variations of the O 2p variance during the ferroelectric phase transitions of 
BaTiC>3 were found to be very small. This surprizing result was explained in terms of the 
electronic structure of this compound as it is interpreted in the Harrison model. 

We think that, when combined with Born effective charges, the band- by-band decomposi- 
tion of the localization tensor could provide a powerful tool for the qualitative characteriza- 
tion of bonds in solids. However, more studies are needed, for different classes of materials, 
in order to make it fully effective. 
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The Born effective charges are in general compared to an isotropic nominal value that is the 
charge expected in a purely ionic compound. All deviations with respect to this reference nominal 
value are referred to as anomalous. 



APPENDIX A: OPTICAL CONDUCTIVITY 

The optical conductivity (imaginary part of the optical dielectric tensor) of a given ma- 
terial is related to its absorption coefficient, the probability of the valence electrons to 
perform optical transitions to the unoccupied conduction bands under the influence of an 
electromagnetic field. If we consider only "vertical" band-to-band transitions (thus neglect- 
ing elementary excitations like the electron-hole interaction or the electron-phonon coupling) 
this quantity writes in the dipolar approximation!^ 

£ a/?M = I 2ft Yl Yl / 7o-^Pnm( k )pL( k )^ (^mn(k) - W) (Al) 
n=l m=N+l JBZ V ; 

where m is the electron mass, p nm (k) = — ih(tp nk \ V^ mk ) and frw mn (k) = £ mk — e nk . The 
matrix elements of the momentum operator can equivalently be expressed as p nm (k) = 
TflUJ nrn (k)(?i„ k |9 k M mk ). 

It has been shown by Souza, Wilkens and Martin0 that e" is related to the localization 
tensor by the relation 

£'U")du=— f ^—{r a r (3 ) c . (A2) 
In order to see the effect of the band by band decomposition, we will write e" as 

= E \ £ 'U^ Gi) + $>'^ Gh Gj) (A3) 

where 

r2„2 



£ 'ap(u;Gi) = 477 6 V V f ^f^p^ m (k)p^ n (k)5 (u mn (k) - u) (A4a) 

el^G^Gj) = 4?r e V V / -^-\p^ m {k)p^ n {k)S (w m „(k) - w) . (A4b) 
The first sum of Eq. ( A4a|) has to be taken over the bands of group G,i while the second 



sum extends over all bands (unoccupied or not) except those of group G^. In Eq. (|A4b| ) , 
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the two sums extend over the bands of group Gi and Gj. It is easy to show that Gj) 
and s" a p{oj; Gi, Gj) are related to the variances and covariances by the relations 

J™ el p {u- GJdco = (r a r,) c (G t ) (A5a) 

J™ e" ap {u- G h Gj)duj = (r a r p ) c (G h Gj). (A5b) 

Thanks to these definitions, the physical meaning of the covariance becomes now obvious: If 
the total localization tensor was simply the sum of the variances (r a rp) c (Gi), the expression 
of the dielectric tensor ( |AT| ) would not only contain transitions between occupied and unoc- 
cupied states, but also transitions between occupied states themselves. It is by adding the 
covariances (r Q r / g) c (Gi, Gj) that one compensates the effect of these forbidden transitions 
in order to get a physically correct quantity 
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